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ABSTRACT 

The steady motion in sunspot magnetic regions is considered 
for both current-free and force-free configurations. The 
sufficient condition for stability is obtained in the presence 
of both external current-free and force-free magnetic fields 
and a steady motion. It is shown that the pattern of such 
steady motion is most important in triggering an instability 
of sunspot magnetic fields, both for the current-free and force- 
free configuration. When there is no steady motion, the current- 
free configuration of sunspot magnetic fields is always stable, 
whereas the stability in the case of force-free magnetic fields 
is connected with the configuration. The onset of a solar flare 
seems to be associated with an instability connected to the 
steady motion within the sunspot magnetic regions. 


*NAS-NRC Associate with NASA 


INTRODUCTION 


It has recently been shown that most solar proton 
flares take place within the &-type sunspot groups (e.g., 

Sawyer , 1965). Furthermore, Sakur ai (1967, 1969) suggests 
that the twisting motion of sunspot magnetic lines of force 
is related to the origin of solar proton flares. He also 
discusses the relation of the 6-type configuration with the 
stability of sunspot magnetic fields. 

The configuration of sunspot magnetic fields has been 
thought to be force-free (e.g., Parker , 1957; G old and Hoyle , 
1960; Longmire , 1963; Gold , 1964, 1968), but there is no valid 
supporting observational evidence. However, it is estimated 
that the magnitude of the electromagnetic force is much greater 
than that of the pressure and gravitational forces within 
sunspot magnetic fields ( Lvist and Schluter , 1955; Gold and 
Hoyle, 1960). Accordingly, the hypothesis that the configuration 
of sunspot magnetic fields is force-free appears reasonable. 

As has been shown by Sakurai (1967), the steady fluid motion 
within sunspot groups, which is connected with the convection 
near the photosphere, must be important in generating an in- 
stability related to the origin of solar flares. 
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In this paper, we will first discuss ste?dy motion 
within sunspot groups and some problems related to the con- 
figuration of sunspot magnetic fields. Then, the stability 
of a force-free configuration of sunspot magnetic fields will 
be discussed. A brief discussion will also be given of the 
mechanism of solar proton flares. 

1. THE STEADY MOTION IN SUNSPOT MAGNETIC FIELDS 
AND THEIR FORCE-FREE CONFIGURATION 


The types of steady motion in sunspot magnetic fields 
will be considered. Fundamental equations appropriate to 
this problem are 


n m (v *v) v = - vP 
c e e,o e,o 


e , o 


ne(E + — v xH)+n m g (1) 
o o c e,o , o o e 


and 


n m. (v. *7) v. *» - 7p + 

o i i,o i,o i,o 


n e(E + — v. x H ) + n m. g (2) 
o o c i,o o o i 


v • n v =0 
o e , o 


7 • n v . - 0 

o 1,0 


(3) 

(4) 


— * 4tt — » — > 

7 x H * — n e(v, - v ), 
o c o i,o e , o 


(5) 
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where v, m, n , E, H, P, g and c are the velocity, mass, 
number density, electric and magnetic fields, pressure, 
gravity force and the speed of light, respectively. The 
subscripts, i, e, and o denote ion, electron and steady state. 
The collision term is neglected because the mean free path 
of particles is so large that the medium is highly conductive. 

The gyro-frequency of ions is denoted as q 
( = eH^/m^c) and the characteristic length of the sunspot 
group is designated as L. If the condition 


| v | , | v | « 0. L 

1 i,o ' e,o i, o 


( 6 ) 


is satisfied, the non-linear terms of eqs. (1) and (2) can 

be neglected (Namikawa and Matsushita , 1970). With the values 
2 9 

of H q - 10 i and L = 10’ cm, the Order of magnitude 

13 -] 

of 0 . L is 8.8 x 10 cm sec which is much larger than c. 

1,0 

Thus, v/e can always neglect the non-linear terms in eqs. (1) 
and (2) . 


Thus, 


-vp 

e,o 


- i - - 

n e(E -f — v x H ) + n m 
o o c e,o o o e 


g 


- 0 


(7) 


-VP + n e(E + - v x H ) + n m g = 0 (8) 

i,o o o c l , o o oi 
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By adding these two equations, we obtain 


VP + — J xH + D g ~ 0, 
o c o o o ’ 


( 9 ) 


where the total pressure, P - P J + P , the mass density, 

o i,o e,o 

p - n (m. + m ) and the electric current density, 
o o i e 

j * n eCv, - v ). The above equation gives the hydrostatic 
J o o i,o e,o 

equilibrium condition. 


When we solve eqs. (7) and (8) with respect to v and 

v? p sJ 

v. , respectively, we have the following equations: 

jL f O 


— ♦ n h «♦ n 1 

v = ” (E x H ) + — =• ~ VP x H 
e,o 2 o o R 2 n Q e e,o o 

o o 


( 10 ) 


and 


cm 

- ~§ « X fi 0 X “ e H 0 
eH 

o 


-4 p — * — • f 1 “ * 

v “ -s (E x H ) - -s — VP x H 
i , o „2 o o „2 n e i,o o 
n rl o 

o O 


cm i-- 

+ —$ g x H o + a i "o’ 


eH 


( 11 ) 


where and a g are arbitrary scalars. Thus, the electric 
current density is given by 
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j « n e(v. - v ) 
°o o i,o e,o 


H 


(7P 


H o ) 


p c ^ 

-§- (g x H ) + n e(a. 

H O O 1 


a ) H 
e o 


( 12 ) 


On the right hand side of this equation, the first two terms 
give the electric current perpendicular to the external 
magnetic field H q , whereas the third term gives the electric 
current along H q and is recognized as a force-free electric 
current. It is clear that this electric current j is in- 
dependent of the electric conductivity of the medium treated 
here (Cowlin g, 1953; Sakurai , 1968). 

If we define the electric currents perpendicular to and 
— « — • — * 

along H q as 3 q Jl and j Q ^ respectively, eq. (5) can be expressed 
as 


v x H 


4tt 

c 


(j O, 1 + ^0,11 ^ * 


(13) 


If the contribution from the pressure and gravitational 
forces is negligible as has been numerically 
estimated (e.g., Ltlst and SchlUter, 1955; Gold, 1968), 
will be negligible. This leads to the conclusion that the 
general steady motion in sunspot magnetic fields is a force- 
free motion: 
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namely, 


v x H * — n e(a, - a ) H . 
o c o i e o 


(14) 


In this case, it follows from eqs. (10) and (11) that 


v ~(E xH)+a H 

e,o 2 o o eo 
n 

o 


(15) 


and 


v. - -s- (E x H ) 4- a, H . 
i,o 2 o o i o 
n 

o 


(16) 


■4 

We furthermore obtain, from the Maxwell equation 7 x E - 0, 

v x (v x H ) - 0, (17) 

o o 

where v q is the mean plasma velocity perpendicular to H q , 

which is easily calculated from eqs. (15) or (16) and equal 
2 -* -* 

to (c/H q ) (E q x H q ) in the present case. The continuity equa- 
tions (3) and (4) reduce to 


v* (n rt 

v + a, 

n H ) = 0 

(18) 

o 

o i 

o o 

v* (n 

o 

v + a 
o e 

n H ) - 0 
o o 

(19) 


The equations obtained above ((14) - (19)) are applicable to 
the steady motion associated with the force-free configuration 
of sunspot magnetic fields. 

If a. * a ( = a) , we have a current-free motion determined 
l e 9 




by the following equations: 

wrr* 

v X H ”0 (20) 

o 

7 X (v x H ) - 0 (17) 

o o 

and 

7(n v -fan H ) - 0 (21) 

o o o o 

In this case, the configuration of sunspot magnetic fields 
is derived from a scalar potential. 

It is likely, however, that the configuration of sun- 
spot magnetic fields is of the force-free type (e.g., A nzer , 

1968; Zwann , 1968). In order to deal with the stability problem 
of both current-free and force-free magnetic fields, we must 
consider the results obtained above for the steady state 
and their perturbation. 

2 . HYDROMAGNETIC STABIL ITY 
2-1 General 

The hydromagnetic stability of the steady motion of 
a plasma will be considered under the influence of an external 
magnetic field. The equations of motions appropriate to this 


problem are 


'IV _ 

|a im$ we* m4 

~r— + (V ,x 7) V + (V * ?) V 

o t e , o e e e,o 


e «. 1 l " ”♦ 

— (o + - v x h + ~ v x H ) , 
m c e,o c e o’ 

e 


*v t _ - 

+ (v. *v) V, + (v •?) V 

dt l,o i i i,o 


—4 *1 ^ 1 —« mm* 

— (e + — v. x h + — v. x H ), 

m^ c i,o cl o 


( 22 ) 

l 


# 


(23) 


where v. , v , h and e are respectively ion and electron 
1 6 

perturbed velocity, magnetic and electric field of the 
perturbation. The second and third terms on the left hand 
side of eqs. (22) and (23) can be neglected in comparison 
with the second and third terms on the right hand side, when 


e, o 



«* L ft 


i,o 


( 6 ) 


|v e l, IvJ « L n t (24) 

and 


We assume here that the change of magnetic field strength 
associated with a solar flare is of the order of li as a lower 
limit . Even in this case, ft - 1.8 x 10 rad/sec and then 
C i L - 1.8 x 10 1 cm/sec. Since both |v e l and | v i | must be 
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smaller than c, the condition (24) is always satisfied. 


Using the Maxwell equation 
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x e * - 


1 

c 


b h 
bt’ 


( 25 ) 


we obtain hydromagnetic equations from eqs. (22), (23), (5) 

•>4 

and (14) where small terms proportional to (b/b t) v x h and 
(v x h) x h are neglected (e.g., Spitzer , 1962; A1 f\6 n and 
Fa It hammer , 1963) as follows: 


-* n e 

Mil 1 — « 

_«_(7xh)xH + (a - a ) H x h 

o dt 4n o c i e o 


(26) 


AU —♦ —4 — < 

•rrr « 7 X (v X H ) + v x (v x h) 

ot o o 


(27) 


or 


.+v x h + v x H =*0, 
o o 


(27) 


where 


r (ra i "i + m e V 
o 


(28) 


The second term of the right hand side of eq. (26) describes 
the effect due to the force-free magnetic field. 

By differentiating eq. (27) with respect to time t, 
and making use of eq. (26), we have 
2 - 

jr . — 7x{[(Vxh)xHJxH + 

. .2 4rr p o o 

ot o 


4tt n e _ ^ 

(c^ - a e )(H o x h) x H o ) + v x (v Q x ^) . 
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( 29 ) 


Let us assume that v and h are proportional to exp (ut) . 
Scalar multiplication of eq. (29) with h and integrating over 
a volume V enclosed by a surface S, gives 


2 f E : 


os I h dV * 
v 


- u> 


f {hx(vxH)}*dS 


-/ 


. r\ 2 r? 2 

(v x h) H q 


dV + uu 


j v Q *{h x (v x h) } dV 


- ID 


/ h x (v x h) • dS 
•i o 


n e f _ 

+ -2- (a - a ) / h-v x {(H x h) x H } dV, 
peter. o o 


(30) 


where we have used the identities 


{ W v * d 


h* v x 


and 


v-{h x (v x h)] « (v x h)*v x h - h-v x (v x h) 
o o o 


together with eq. (26) and Gauss's theorem. According to eq. 
(27) ’ ,the first and forth terms on the right hand side of eq. 
(30) are reduced to 
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uj / (h x e)*dS 
J S 

In the outer space W - V of the sunspot magnetic regions, 
where W denotes the whole space, we obtain 


( 31 ) 


/ v.(h x e)*dV * - / (h x e)-dS 

S 


-/ 


W-V 


f e*v x h dV - / h*Vx 

w-v w-v 


e dV 


/ 


a) / h dV 
W-V 


(32) 


where we have used eq, (25) and the condition 


Vxh-OinW-V 

This indicates that perturbation current flows only in the 
volume V. 

Eq (30) is thus reduced to 

h ^ + h “ + X 3 = °’ 


where 




dV 


{ (Vx h) x h} dV 


(33) 


(34) 


(35) 


(36) 


and 
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* 


A 


V 4n P o 


t 2 J2 

<7 x h) *H dV + 
o 


f — (a. - a ) H 2 h-7 x h dV, 
•i. c c i e o 


V p o 


(37) 


where we have used the identity 


h-Vx[(H x h) x H } “ H h-v x h 
o o o 


(38) 


The first term of I g is always positive, but the sign of 

aw* 

the second term is dependent on (a. - a ) h*v x h. 

i e 

In order to study the necessary and sufficient condition 
for stability, we solve eq. (34) with respect to w : 

arr^iV 1 '- 41 ! V 


0) 


(39) 


Thus, the necessary condition for stability is 


i 

, t 

'■u 

t. * 


I 2 > 0 and - 4 I 3 ^ 0 


or 


(40) 


I 2 » 0 and I 3 > 0. 

Under the condition (41), we have a pure oscillation of 


(41) 


frequency <« given by 


U> 


41 


(42) 


. -' Jl 


2-2 Current-free fields 


In this case, the second term of eq. (37) is neglected 


since 


V x H = 0 
o u • 


( 20 ) 
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It thus follows that 1^ is always positive. Consequently, 
the condition for stability is g;iven by eqs. (40) and (41). 
In order that the current-free magnetic configuration is un- 
stable, the following condition must be satisfied: 


I 


2 



{ (v x h) x h} dV < 0. 


( 43 ) 


If we use a dyadic notation (J ackson , 1962), this result can 
be written as 


where 


and 



(7-T) dV, 


h h - 


1 h' 



, e + e„e_ 
11 2 2 


+ W 


(44) 

(45) 


e , e and e denote the unit vectors of Cartesian coordinate. 

X. J 

In general, the equation of momentum change in a magnetic field 
is given by 


— * 

dP 

dt 


7-T 


Therefore, the term v Q (dP/dt) (-v^- (7 -T) )expresses the rate of 

work due to the motion v . Thus, I n gives the rate of the work 

o 

done by the steady motion v q for the whole space where this 
motion is observed. 
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The inequality (43) shows that the above rate must be 
negative. This means that the initial magnetic state is 
not at energy minimum, but tends to move to a lower energy 
state whenever a perturbation is introduced. We can estimate 
the sign of the integral Ig by considering both modes of magnetic 
perturbation and steady motion. We can conclude that the pattern 
of steady motion is most important in triggering an instability 
of the sunspot magnetic fields. 

2-3 . Force-free Fields 

When eq. (37) is not negative, the criterion for 
the stability of current free magnetic fields can be applied 
to that for force-free magnetic fields. In this case, when 
I 2 < 0, the configuration of force-free fields is also unstable 
provided there is some steady motion. 

When there is no steady motion (and therefore 
I 2 =0), we obtain another criterion for the instability of 
force-free magnetic fields such as 

I 3 < 0. (46) 

Consequently, the following inequality must be satisfied: 

/ 2^ 2 r n e 2^ 

7 (7 x h)‘H d V + / - 2 - (a. - a ) H h- (7 x h) dV < 0. 

4n p o J p c l e o 

vo Vo 

(47) 

Since the first term is always positive in this equation, the 

second term must be negative and its absolute value must be 
greater than the first term. 
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3. IMPORTANCE OF STEADY MOTION ON THE STABILITY OF 

SUNSPOT MAGNETIC FIELDS 

As discussed in the last section, the steady fluid motion 
plays an important part in generating the instability of sun- 
spot magnetic fields, both current-free and force-free. As 
we know, the growth and decay cf sunspot magnetic fields is 
well connected with the convective flow near the photo- 
spheric surface (e.g., Parker , 1955; Babcock , 1961; Sakura i , 
1967). The pattern of this flow seems to be most effective 
in generating the instability of sunspot magnetic fields. 

At present, the configuration of sunspot magnetic fields 

is thought to be force-free (e.g., Lust and Schluter , 1955; 

Gold and Hoyle , 1960; Sakurai , 1967; 1969; Anze r , 1968). 

If the induced force-free current |n e(a. - a ) h| is much 

o i e 

smaller than the induction current |c (v x h)/4tt|, the effect 
of the ambient force-free current seems to be negligible in 
the problem of the stability of sunspot magnetic fields 
even if their configuration is force-free. Therefore, the 
criterion for the stability given by eqs (40) and (41) can be 
applied to the problem of the stability of force-free magnetic 
fields. The fluid motion is both necessary and sufficient 
to produce the instability of both current-free and force-free 
magnetic fields as discussed. 
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The development of sunspot magnetic fields is associated 
with the convective motion near the photosphere (e.g., Parker , 
1955; Babcock , 1961; Sakural , 1967). This motion may be 
responsible for the generation of sunspot magnetic instability. 

An instability produced in this way seems to be connected 
with the triggering of a solar flare. As suggested by Sakurai , 
(1967, 1970), the ambient convective fluid motion may be 
associated with the rotating motion of sunspot magnetic regions 
and the type of sunspot groups. The fact that the onset of 
solar proton flares is well correlated with this rotating 
motion suggests that some steady ambient fluid motion is 
very important in triggering solar flares. 

4. CONCLUDING REMARKS 

It has been shown that the steady motion plays an important 
part in the stability of sunspot magnetic fields. Such motion 
may be associated with the growth and decay of sunspot groups 
and also related to the convective flow near the photospheric 
surface within sunspot magnetic regions. 

The steady motion within sunspot magnetic regions and its 
pattern must be studied in detail when we deal with the stability 
problems of sunspot magnetic fields. The results obtained here 
may have application to other astrophysical problems. 
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